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1 .  Introduction : 

Stokes  (1846)  showed  that  a  sloping  beach  can  act  as  a  wave 
guide  which  allows  gravity  waves  to  travel  along  shore  with  an  ampli¬ 
tude  decreasing  exponentially  with  a  distance  from  the  shore.  Since 
then,  several  authors  including  Ursell  (1952),  Greenspan  (1956)  and 
Hunk  and  others  (1956)  discussed  long  or  short  gravity  waves  on  a 
sloping  beach  of  a  homogeneous  water  either  from  a  mathematical  point 
of  view  or  from  application  to  geophysical  problems. 

Kelvin  (1879)  treated  a  problem  of  long  waves  propagating 
along  a  straight  coast  in  a  rotating  sea  and  found  that  the  amplitude 
decreases  exponentially  from  the  coast  to  the  left  hand  side  of  the 
direction  of  propagation  in  the  northern  hemisphere.  Energy  of  such 
waves  which  are  called  Kelvin  waves,  therefore,  is  seemingly  trapped 
by  the  coast. 

The  edge  waves  over  a  sloping  beach  in  a  rotating  sea  of 
homogeneous  water  was  treated  by  Reid  (1958)  and  Kajiura  (1958).  They 
showed  that  in  thi3  system  there  are  two  inert io- gravitational  waves 
travelling  in  opposite  directions  and  one  quasi-geostrophic  wave. 

However,  density  stratification  of  the  water  as  well  as  a 
beach  slope  becomes  important  for  problems  concerning  response  of  the 
marginal  sea  to  atmospheric  disturbances  and  for  problems  concerning 
internal  waves  near  the  coast.  In  order  to  represent  such  stratifica¬ 
tion  in  a  way  which  is  feasible  to  mathematical  treatment  and  yet, 
keeps  essential  dynamical  feature  of  the  system,  a  two-layer  model  of 


2 


Che  thermoclina,  consisting  of  two  fluids  of  different  densities  Is 
adopted.  Deviation  of  modes  of  wave  motion  In  a  two-layer  model  from 
those  In  a  continuous  density  was  discussed  by  Eckert  (1960)  In  a  more 
or  less  general  way. 

Essential  features  of  motion  near  the  coast  with  a  length 
scale  of  the  width  of  a  continental  shelf  and  with  time  scales  longer 
than  a  few  hours  can  be  represented  by  a  mathematical  model  for  waves 
In  a  two-layered,  rotating  sea  bounded  by  a  straight  coast  and  with 
a  variable  depth.  There  are  two  types  of  bottom  topography  which 
renders  mathematical  analysis  feasible  and  yet  simulates  an  area  close 
to  the  continent. 

One  type  is  a  continental  shelf  with  a  uniform  or  variable 
depth  connecting  with  the  open  sea  with  a  uniform  depth  as  treated  by 
Ichiye  (1963).  The  other  type  of  simple  bottom  topography  is  a  depth 
increasing  linearly  from  the  coast.  In  a  two- layered  sea,  there  must 
be  longshore  geostrophic  currents  due  to  pressure  gradients  if  the 
Interface  as  well  as  the  surface  intersects  the  bottom  at  the  coast. 

This  study  has  treated  such  a  model  mathematically.  In 
order  to  simplify  the  analysis,  it  is  assumed  that  geostrophic  currents 
in  the  upper  and  lower  layers  have  uniform  velocities  owing  to  (he 
surface  elevation  increasing  linearly  from  the  coast.  Although  such 
an  infinite  uprise  of  the  sea  level  seems  to  be  unrealistic,  the 
waves  considered  have  amplitudes  significant  only  near  the  coast  and 
the  simplified  model  may  be  justified  if  Interpretation  of  mathema¬ 
tical  results  is  proper.  Slopes  of  the  surface  and  the  Interface  are 
assumed  so  small  that  the  approximation  of  hydrostatic  pressure  is 
valid  as  shown  by  Eckart  (1951). 
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2.  Fundamental  Equations : 

The  x-  and  y-  axes  are  taken  on  the  level  surface  perpendi¬ 
cular  and  parallel  to  the  coast-line,  respectively,  and  the  s-axis 
is  taken  vertically  downwards,  as  shown  in  Fig.  1. 


Fig.  1  Schemetic  diagrams  of  coordinate  system  (A)  and  a 
vertical  section  (B). 

It  is  assumed  that  densities  of  water  in  upper  and  lower 
lr.yers  are  uniform  and  their  difference  Is  small.  The  basic  currents 
ar<e  geostrophic  and  their  velocities  are  constant  in  each  layer.  The 
waves  are  considered  to  be  perturbations  to  such  basic  currents.  The 
non-linear  inertia  terms  are  neglected  with  assumptions  that  the 
vertical  velocity  and  amplitudes  of  the  waves  are  very  small.  The 
vertically  integrated,  linearized  equations  of  motion  without  dis¬ 
turbing  forces  and  equations  of  continuity  are  respectively  as 
follows: 
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In  these  equations  suffices  1  and  2  correspond  to  the 
quantities  in  the  upper  and  lower  layer,  respectively;  Mj,  Nj  are 
horizontal  mass  transports  respectively  along  the  x-  and  y-  axes, 

1 7j  i  the  basic  current  assumed  to  be  geos  trophic ,  f  ^  the  density, 
hj  (*SjX)  the  thickness  of  the  layer;  ^  ^  and  4>  2  are  respectively 
the  depression  of  the  surface  and  tha  interface  from  the  equilibrium 
levels,  ("  y  2"  f  density  difference,  £  the  Coriolis  parameter, 
£L  the  gravitational  acceleration,  and  ^7  horizontal  gradient  operator. 

The  assumption  of  geostrophic  balance  for  the  basic  current 

leads  to 


Assuming  that  horizontal  mass  transports  and  displacements 
of  the  surface  and  the  interface  are  proportional  to  exp£i(&|J  -<T £)} 
the  equations  of  motions  (1)  becomes  a  system  of  linear  equations  for 


and  N ,  with  constant  coefficients. 

In  tr:  .as  cf 


are  given  by: 


The  solutions  of  Mj  and  Nj 
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Substituting  Mj  and  Nj  into  the  equation  of  continuity  (2),  we  can 
get  equations  for  ^ (  and  *• 


l»-(w  *$)!<■ 


(6) 
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and  , 

(7) 

where 

Aj  =  .  Sj  -  dfcj/dz 

and 

&  ~  ^  Z]  ■ 

Elimination  of  ^rom  aquations  (6)  and  (7)  yields  the  equation 

for  : 

([>-  a.  ID  +b)$,  =0  , 

(8) 

where 
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The  equation  (8)  can  be  separated 
Laguerre  differential  equations: 

(D-  K,)$>  = 0 

or  CO  -  Kz)  =  0  , 

"here  K,  =  Ha  ±  I^U  )  .  <l3> 

It  Is  seen  that  a^  -  4b  ^0  for  small  S'  . 

In  order  that  the  solutions  of  equations  (11)  and  (12) 
may  be  finite  both  at  the  coast  and  at  the  Infinite  distance,  the  Kj 
must  equal  to  *(n+  l)k,  where  n  Is  an  integer.  These  conditions 
yield  the  relation  between  £  ^and  £"2  and  the  frequency  equation, 
respectively: 


Into  two  independent 
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(12) 
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(14) 


and  (*>”  k*V,Wi  4  (2**0 


2(‘0,+Ia)z')  + 

+  f^WiWiO-cuj)-r  &}u/jU/,0-wS)}J  +  $&l*£z+ 

+  •“7 £oua(/-ivl)+  $&&u/t(l-io,*)  4 60,C*J2(i-4S,9Xt-tOf)*0 


(15a) 


The  frequency  equation  (15a)  can  be  reduced  into 

tlj  +U,(I  -U?  )  J  f  (tv+Mz  4 1 }  4  w*(  |-  <oa*>  3 


=  0-6  '>0j,(D3O'-to'*X\-U)f') 


(15b) 


In  these  equations  following  notations  are  used: 


and 


OJj  S=  5u  *  CU-  sC^'l 


7§*» 


/*  i*  1  l 

The  equation  (15b)  can  be  written  In  the  polynomial  about  w.  The 

coefcicients  of  different  powers  of  w  are  shown  in  the  following 

Table  1,  In  which  the  dimensionless  quantities ,  £  f£,  and  Y are 

defined  by  -  \£  «/  VT& 

1  = -  i  - k  -  — 

Equation  (3)  yield.  (|+£)fe,  =  ) Vf 1  =  jSCS.tSOM'1 


v-  _  Si  4  Sa 
*  ~  “ 


Therefore  (I4£)  ^  Is  always  positive,  no  matter  in  what 
direction  the  current  In  upper  and  lower  layers  flow.  Values  of 
fox  different  values  of  and  wave  length  are  shown  in  Table  II, 
fox  -f  ■  l(f4  (sec”^). 
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Table  1  Frequency  Equations 

(J1 _ CoefKcierits _ 

U)4  / 

cu5  -30-£)&' 

cos  i2)i'st (wnxttvfftt-irxf-Dta-o \ 

rlso-t)  -(iTtyffi  fc' 

CO2  -3(f~3Z  ti2)  tt?  f  *(*9»  r-l  -if  +i*  )fc3+ 

+  Ot€)j-f  +  3(r-l -£)}&'*+  '+  I 


CO '  - J ix( /-£)&' r-  (J*i-<X»£)jrf(Y-iX3E- 1)  +El(£-3)}K*- 

-(hi'>jrl(2*Ttx/-t*xv-\)  +3tt*+r-i)}  &%x*m i/ftsrtf-Ofi  A'*+ 

+  Oti'fj'Tk' 

+(wHXftvf?r/t£')U* 

-Cr*tbx(tTif£(r~i)±-£/4i-  (2*+>X/+£ )j-Uf-£'XH )r£T}&  ^ 
+-£(n*)fcTfc'* 

o 

Table  II  Nondimensional  wave  number  "fL  —  Vfe  / *f  corresponding 

to  a  given  wave  length  L  (km)  and  velocity  V,  (m/sec) , 
where  f  ■  10“^sec"l. 
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3.  The  waves  for  n  -  0 

When  n  -  0,  the  polynomial  (15b)  can  be  factorised  as  follows 

&/j-  0  —  (*+£)&'  }i(tu+£fk.'Xtti1'£&“  •)  *•* X 8 
+  Ci-  syco-  &Xcv+£(tiXu>-{L'+iXiv+ffi'-*j>j(i6) 

It  is  seen  at  once  that  equation  (16)  has  two  roots: 

60=  *'-|  (17) 

CO  =  -£&'-{  (18) 

The  root  (17)  or  (18)  makes  the  denominator  f  — <Jy  in  equations  (4) 

amd  (5)  vanish  for  j  -  1  or  2,  respectively.  Therefore,  equations 

(.4)  and  (5)  yield  m  „-AX 

(j  ^  e 

in  which  the  suffix  j  equals  to  1  or  2,  corresponding  to  the  root 
(17)  or  (18),  respectively. 

Substituting  Yj  Into  cither  (6)  or  (7),  it  is  found  that 
and  are  proportional  to  £ "  .  The  mass  transport 

or  M2  respectively  for  the  root  (17)  or  (18)  can  be  derived  from 
one  of  equations  of  motion  (1)  and  equations  of  continuity  (2). 

The  boundary  condition  that  Mj  thus  determined  should  vanish  at 
the  coast  and  be  finite  at  an  Infinite  distance  leads  to  a  relation 
which  determines  the  values  of  k. 

One  such  corresponding  to  the  root  (17)  is  given  by 

+2tyiaj,k  +2f+W,ttix0-cu*)W'> 

Since  w^  ■  -1  and  W2  is  linear  with  k,  equation  (19)  is  a  cubic 
equation  about  k.  This  equation  yields  the  roots  when  £  is  very 
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mail. 


and 


i 


+  Oti)  , 


A  — -  +  oo) , 

6(3S,+b)3 


4 


&<%f§ 


Relation  (20a)  Is  equal  to  2f  */%($, +S* ) 


(20a) 

(20b) 

(20c) 


when 


;«  0  and  •  Vg  ■  0,  as  obtained  by  Reid  (1958). 

The  waves  corresponds  to  tie  roots  (17)  and  (18)  are  degen¬ 
erate  cases  of  other  modes  of  waves  as  In  case  of  a  uniform  density. 
Since  these  waves  are  possible  only  for  discrete  wave  mxnbers ,  they 
cannot  be  excited  by  local  atmospheric  disturbances  of  arbitrary 
dimensions . 

Other  four  roots  of  w  in  equation  (16)  are  determined  from 
the  following  biquadratic  equation  with  w: 

CV*  -  2i«-*)k’+i}  w3-  £ *£K'\(H*MM/-36-e)A-if<0* 


-*  \-2(z- -  a  -  At +»)£'  V  a(tt£Xr-i-s)&'tA'-o-z) l'+ 
*  oKiim' }  w  +  Ce*&^£(s-ok'3-(iT2x*-  i+zWi'3- 

-  £&*  *  o 


(22) 

As  seen  from  Table  II,  "£' «  $  "l  for  the  wave  lengths  longer 
than  several  kilometers  and  reasonable  values  of  because 
In  such  a  range  of  ft  roots  of  equation  (22)  can  be  expressed 
epproxlmately  with  algebraic  relations  of  A 

First,  only  terms  with  S'1  are  kept  for  the  same 
powers  of  fof  ,  since  %  «l  .  Further,  if  £  <</  ,  powers  of 

higher  than  the  first  in  the  coefficients  of  jg2  and  w  can  be 
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neglected.  The  left  hand  side  of  equation  (22)  can  be  written  as 
Cu(£u-i)£co2-1  Co  -  (/tOS^fc'] 

(23) 

-f  (t+ixt-nyLg3 *o 

If  the  last  two  terms  of  equation  (23)  are  neglected,  the  first 
approximation  of  four  roots  of  (23)  is  given  by: 


CO  O  CO  /  ,  (24a)  (24b) 

and  ,  ,,  // 

<u * ±  *  i i  +  ft'}  ,  <25i)  <25b) 

In  relations  (25a)  and  (25b),  %  -terms  without  a  factor  S'1  are 
neglected.  The  following  second  approximations  for  (24a)  and  (24b) 
arc  determined  from  (22)  to  the  order  of  fe2  and  -%_/ 1  respectively. 

^  -fr+  Lining 

0 

Cv  -v  /  -f  l  I  +  Cl  —  2. IT  1  & 


(26a) 

(26b) 


The  term  of  -fa1  in  (26a)  becomes  important  when  2  =  0  .  The 

roots  (26a)  and  (26b)  yield  respectively  following  phase  velocities: 


C  - 


(27a) 


C  =  77*  (f  +  ^toR)  l  ) 

The  depression  of  the  interface  corresponding  to  (26a) 
(26b)  are  respectively  given  by 

^  '  ’f  T^2 


(28a) 

(28b) 
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Sine*  aquations  (28a)  and  (28b)  Indicate  that  displacements  of  the 
Interface  are  much  larger  than  and  reserve  in  direction  to  those  of 
the  free  surface,  the  waves  corresponding  to  roots  (26a)  and  (26b)  are 
Interpreted  as  the  barocllnlc  mode.  The  frequency  corresponding  to 
(26a)  Is  very  small  and  the  mass  transports  which  are  obtained  by 
substituting  (26a)  into  (4)  and  (C)  are  quaslgeostrophlc  or  geostrophic 
to  the  order  of  magnitude  of  ik,  .  Since  the  frequency  for  the  root 
(26b)  almost  equals  to  that  of  Inertial  oscillations,  waves  corres¬ 
ponding  to  this  root  may  be  called  inertlogravitational  waves.  The 
phase  velocity  (27a)  of  quaslgeostrophlc  waves  shows  that  the  waves 
are  carried  by  the  current  of  the  lower  layer.  When  0,  they 
propagate  to  the  positive  y-directlon.  When  V^?0#  they  may  propagate 
to  the  negative  y-  direct  Ion  for  smaller  but  to  the  positive 

y- direct  ion  for  larger  /£/  in  a  range  of  satisfying  £<</  •  The 

value  of  at  which  such  change  In  propagation  direction  occurs  is 
given  by  £  f(f-t  which  lies  In  the  range  of  H'  considered  if  £ 
is  sufficiently  small.  Inertlogravitational  waves  are  carried  by 
the  current  of  the  upper  layer,  as  seen  from  equation  (27a). 

Dynamical  properties  of  the  waves  corresponding  to  (25) 
may  be  understood  when  we  consider  a  range  of  ;£  satisfying  $.'<<■  Sf*?)  t. 

,  Under  this  condition,  the  roots  (25)  can  be 

expressed  by 


Oy  ^  ~  (It  Ik 

(29a) 

Co  %  / •+  ( H  £)y k 

(29b) 

The  displacement  of  the  interface  corresponding  to  (29a)  and  (29b) 
is  given  by: 

£  _  Jti-  k 
S*  —  Y 


(30) 
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which  Indicates  that  the  waves  are  of  barotropic  node.  The  same 
argument  as  for  the  waves  corresponding  to  (26a)  and  (26b)  can  be 
applied  to  the  waves  of  barotropic  mode  and  waves  for  (29a)  or  (29b) 
may  be  defined  as  quasl-geostrophic  or  inertiogravitational  waves, 
respectively. 

When  &  1 1  +£  I H  1  A  <<  X&  1^1 ,  equation  (22) 

can  be  approximately  separated  into  two  quadratic  equations: 

4)2-2[l+0-S)f/W  -0+oIfc'aO  (31) 

o 

fcu-f2(r-i-2)fe./+> (32) 

The  roots  of  (32)  represent  the  waves  with  high  frequencies  which 
bocome  waves  of  barotropic  mode  corresponding  to  the  roots  (29a)  and 
(29b)  as  ~jk.  decreases  to  the  order  less  than  £  .  The  waves  corres¬ 
ponding  to  the  root 8  of  equation  (3  0  have  lower  frequencies  than 
those  for  (32).  They  become  waves  of  barocllnic  mode  corresponding 
to  the  roots  (28a)  and  (28b)  as  £  decreases. 

The  four  roots  of  (22)  for  a  range  of  "A.  from  10~^  to  10 
are  numerically  computed  by  taking  following  constants. 

1  =  0  9  OS  ,  /  *«<*'*>  >  (33) 

The  curves  of  w  are  plotted  against  in  fig.  2a,  b.  The  values  O 
and  oo  of  £  corresponds  to  the  cases  of  no  current  in  the  lower  and 
upper  layer,  respectively.  Since  vanishes  for^***0,  is  replaced 
by  fa-"  (=  fk'Vif'1 )  in  ^8*  2b.  It  is  seen  that  general  features  of 
four  curves  of  w  are  very  similar  for  different  values  of  £  except 
those  of  quasigeostxophic  waves  of  barocllnic  mode.  Since  the  prom¬ 
inent  terms  in  most  approximate  formulas  for  w  include  (i+Z)%,/ 
as  a  combination,  the  curves  of  w  are  almost  the  same  for  negative 
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values  of  (hi )  as  for  positive  values  f  /  Is  taken  as  an  argument . 

In  fig.  3  the  ratio  j  $z  ( |  for  four  different  roots  of 
the  biquadratic  equation  In  (16)  ere  plotted  against  for  a  case 
of  £  *  0  .  These  ratios  are  almost  unchanged  In  a  range  of 
considered.  Although  explicit  formulas  of  Indicating  aeparatlon 

of  barotroplc  and  barocllnic  modes  such  as  (28)  and  (30)  are  derived 
only  for  a  certain  range  of  H'  ,  such  separation  la  almost  complete 
In  the  whole  range  of  H',  There  seems  to  be  no  important  Interaction 
between  barotroplc  and  barocllnic  modes. 

Since  the  waves  of  barotroplc  mode  are  similar  to  those 
treated  by  Reid  (1958),  the  approximate  formulas  of  the  roots  of 
this  mode  can  be  compared  with  his  results.  From  the  definition  of 
t  ,  /  ,  V*  and  $  ,  we  have 

(»£)!(  k'  «•  &  (34) 

If  we  put  s  =  S'/  +  S'*  ,the  roots  expressed  by  (25)  become  the  same 
as  those  of  equations  (35)  and  (36)  in  Reid's  (1958)  paper.  (Mote 
that  his  notation  w  is  equivalent  to  -<f  or -fu/  In  our  notation.) 

As  jz  or  \  increases,  these  roots  can  be  approximately  expressed  by 
roots  of  equation  (31),  which  are  equivalent  to  those  of  equation 
(40)  In  Reid's  paper. 

Since  (/+£)&'  is  positive  from  relation  (34)  and 
from  the  definition,  the  approximate  formulas  (23),  (31)  and  (33) 
always  yield  real  values  of  g.  This  situation  Indicates  that  the 
waves  for  n  ■  0  are  stable  for  a  range  of  wave  numbers  satisfying 
the  condition  that  o  <  &  « <5H-f  PT_I  .  This  range  Includes  waves 
whose  wave  lengths  range  from  thousands  to  several  kilometers,  for 
ordinary  magnitudes  of  £  (~  10"3)  and  (  ~  //m/s).  Such  per¬ 
sistence  of  stability  in  the  range  of  wave  lengths  is  different  from 
occurrence  of  long  unstable  waves  which  are  predicted  in  a  mathema- 


t 


versus 
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tlcal  model  of  cyclones  es  waves  on  an  Inclined  boundary  of  two  air 
masses  in  the  atmosphere  (Godske  and  others ,  1957).  This  difference 
might  certainly  be  due  to  the  kinematic  constraints  imposed  by  the 
bottom  configuration!  the  free  surface!  and  the  coast  In  the  present 
model,  though  a  real  reason  must  be  yet  to  be  found  In  more  elaborate 
study.  (See  Appendix  A) 


i 

I 

< 


16 


4.  The  waves  for  ft  2  1 

Whan  A,'  la  much  smaller  than  S(27t+f)m/  ,  equation 
(15b)  yields  approximate  relations  for  CJ  •  Approximation  for 
three  roots  in  a  low  frequency  range  is  obtained  by  equating  to  sero 
respectively  the  linear  and  quadratic  part  about  co  in  equation 
(15b),  in  which  only  the  seroth,  first  and  second  powers  of  A. /  *re 
retained  respectively  in  the  coefficient  of  the  second,  first  and 
seroth  power  of  (O  in  the  polynomial  shown  in  Table  1.  These  roots 
are: 


co — (i-£*x»-or,+E$&4  <35*) 

(J'ill  TH  (35b) 

As  the  roots  (26a)  and  (26)  for  71  •  0,  these  roots 
represent  motion  of  barocllnlc  mode  which  is  characterised  by  the 
valuas  of  ^  proportional  to  £,  .  Substitution  of  (35a)  and  (35b) 
into  (4)  and  (5)  as  well  as  the  magnitude  of  characteristic  fre¬ 
quencies  confirm  that  the  roots  (35a)  and  (35b)  correspond  res¬ 
pectively  to  quaslgeostrophlc  (  JL  ■  6)  and  inert io- gr avlt at ional 
waves  (  £  ■  4  and  5) ,  where  £  denotes  the  mode  of  waves  for  different 
roots  of  (15b)  with  £  •  1  to  3  and  4  to  6  corresponding  to  barotropic 
and  barocllnlc  mode,  respectively. 

The  three  roots  of  (15b)  in  a  high  frequency  range  is 
determined  from  a  cubic  equation: 

-  (  ) 

If  the  term  of  a)zis  neglected,  this  equation  becomes  the  same  as 
equation  (30)  in  the  paper  of  Reid  (1958)  by  putting  S  *  S,  -t  Sx  • 
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Therefore,  these  roots  correspond  to  the  motion  of  barotropic  mode. 
The  approximate  formulas  for  three  roots  of  (36)  con  be  written  as: 

Cj - S-'(h-Z-)  ?&'  (37.) 

co~  1 1  <37b) 

with  the  condition  fe'«  §“  .  The  same  argument  as  applied  to  the 

roots  (35a)  and  (35b)  indicates  that  the  roots  (37a)  and  (37b)  cor¬ 
respond  to  quaslgeostrophic  (  «  3)  and  lnertiogravltatlonal  waves 

{  JLm  1  and  2),  respectively. 

In  a  range  of  satisfying  the  condition  that 
!)-*«  £  /«  &~*(2?1+t )  equation  (15b)  can  be  separated 

into  a  quadratic  and  biquadratic  equation  of  oO  .  The  quadratic 
equation  is 


CaJ2-  (271+1  Yh-£) Y *k'  ~  0  •  <38a> 

If  5  is  svibstituted  for  5,+S*.  and  ^  is  used  instead  of  A.  ,  this 
equation  becomes: 

0Ja —  (29H0 —  O  (38b) 

Therefore  the  roots  of  (38b)  represent  the  gravitational  waves  of 
barotropic  ( l  -1  and  )  similar  to  Stokes '  edge  waves . 

The  biquadratic  equation  of  U>  is  obtained  by  equating 
to  zero  the  sum  of  all  the  terms  with  a  factor  in  Table  1.  If 
g(2 71+i')~l  «  <£  (271+1  f  this  equation  has  the 

following  four  roots: 

Co  ^  (2+1+i  )  1  -  £  ""  (39a) 

to  i  I  •¥  0  ( ) 


Cu  a,  - 


27I+  I 


+  OC*')  . 


(39b) 

(39c) 
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The  roots  given  by  (39a)  and  (39b)  correspond  respectively 
to  (35a)  and  (35b)  and  represent  quaslgeostrophlc  (  £  -  6)  and  lner- 
tlogravitatlonal  -4  and  5)  waves  of  barocllnlc  mode.  The  root 
(39c)  represents  quaslgeostrophlc  waves  of  barotroplc  mode  ( t  ■  3) 
for  -k  >?  f  r*  (p)"1  In  Reid's  (1958)  paper. 

As  £'  approaches  to  and  exceed  (&■*+/)  ~f  ,  roots  of  the 
biquadratic  equation  are  no  more  adequately  expressed  by  approximate 
formulas  (39a)  to  (39c).  The  formulas  for  such  roots  cannot  be 
determined  as  simple  algebraic  relations  of  -fcf  . 

In  figure  4,  the  curves  of  six  roots  of  equation  (15b)  with 
n  «■  1  are  plotted  for  fc/  .  A  range  of  and  constants  &  and  T 
are  the  same  as  for  the  curves  of  figure  2  of  n  ■  0.  The  curves  of 
roots  of  barotroplc  mode  have  features  similar  to  those  of  Reid's 
(1.958)  paper,  Indicating  that  there  Is  little  Interaction  between 
barotroplc  and  barocllnlc  modes  In  a  range  of  jkf  considered.  Also, 
numerical  calculation  as  well  as  approximate  formulas  for  cO  indicate 
that  the  roots  of  equation  (15b)  are  real  for  /n  Z  f  and  the  wa/es 
arc  stable  in  this  range  of  4/-  • 
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5.  Group  velocity. 

The  group  velocity  Cr*  *  (ft)  of  the  waves  corres¬ 
ponding  to  the  root  ajpn  (-ft!)  of  equation  (15b)  is  given  by 


f  /j>/%  dGtn  _  &&!* rr 


(40) 


in  which  GJn  is  the  frequency  corresponding  to  UJf.rt  and  the 
suffices  £  (  -  1  to  6)  and  n  ("0,1,2...)  are  the  mode  end  the  order 
of  the  waves.  If  F (cu,  £/)  represents  the  frequency  equation  whose 
coefficients  of  powers  of  LO  are  listed  in  Table  1,  the  group 
velocity  (40)  is  expressed  by: 


(*)-*{-  (Ij )/(!£) <41) 

General  features  of  group  velocity  of  waves  of  different 
inodes  may  be  seen  from  the  slope  of  the  curve  of  OOg^  (jO  plotted 
against  as  well  as  approximate  formulas  for  Qjfm  (-$_') 

Trends  of  roots  (-fe.')  of  barotropic  mode  as  functions  of 

are  almost  similar  to  those  discussed  by  Reid  (1958)  in  a 
range  of  considered  in  this  study.  Thus,  the  group  velocities 

of  the  barotropic  modes  have  similar  features  to  those  of  his  theory. 

In  order  to  discuss  the  group  velocities  of  barocllnic 
mode  for  n  -  0,  approximate  formulas  of  c O  are  determined  from 
eq'uation  (32)  for  two  ranges  of  .  In  a  range  {>?£/>?$  » 

roots  of  (32)  are  given  by 

Oo  ^  -efe'+fO-r'X n*£) -£**'*}  £*  (42a) 


UJ  ^  I  +  2  ( I- Y”*- Y''?)  &/ 


* 


(42b) 
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In  a  range 


S‘,»n,»i 


,  these  roots  ere  expreeeed  by: 


ooz.  2  1  +  (43a)  (43b) 


where  +0  + 1\\- r')}r» 

%>  *  (i-r-e-ur'Xi^r1  . 

The  formula  with  minus  or  plus  sign  before  the  double  sign  is  denoted 
by  (43a)  or  (43b),  respectively.  Approximate  formulas  of  cO  in  the 
quasigeos trophic  waves  of  baroclinic  mode  are  expressed  by  (26a)  (42a) 
and  (43a)  in  respective  ranges  of  A'  .  Those  of  inertiogravita- 
tional  wave  are  given  by  (26a) ,  (42b)  and  (43b) . 

The  group  velocity  as  well  as  the  phase  velocity  for  the 
quasigeostrophic  waves  is  smaller  than  corresponding  velocities  of 
the  inertiogravitatlonal  waves  in  a  range  of  A '  considered  as  seen 
from  relations  (26),  (42)  and  (43).  If  £  is  finite  and  of  order  of 
magnitudes  less  than  unity,  the  group  velocity  of  the  quasigeostrophic 
waves  is  almost  the  same  as  the  phase  velocity  tfiich  equals  nearly  to 
,  velocity  of  the  permanent  flow  in  the  lower  layer  in  a 
range  «  |  .  Therefore,  with  such  a  condition  these 

waves  are  almost  non-dispersive.  If  Z  ^  0  ,  in  a  range 

H ' «  J  the  group  velocity  of  quasigeostrophic  waves  are  twice 

the  phase  velocity  which  is  the  order  of  magnitude  of  A'  Vi 
In  a  range  S ”*  »  |  ,  group  velocity  is  almost  the  same  as  the 

phase  velocity  for  any  values  of  £ 

The  group  velocity  of  the  inert iogravity  waves  equals 
approximately  to  ZV[  in  a  range  $.  <<  1  »  and 

to  2nTI  V(  in  «  range  / «  &'«  ,  if  Jf  »  |  and 

£  «  j  .  It  is  much  smaller  in  a  range  <<  |  than  the 

phase  velocity  which  is  nearly  equal  to  V' /£'(-  -fA'*),  Therefore, 
these  waves  are  strongly  dispersive  in  such  a  range  of  A'  • 
They  become  less  di  sperslve  in  a  range  |  «.  £ ' «  &  as  in  the 

quasigeos trophic  waves. 
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Analogous  discussions  may  be  applied  to  the  group  velocity 
of  the  waves  of  order  higher  than  seroth  order.  The  group  velocity 
of  the  second  inertiogravlty  waves  corresponding  to  the  root  with  the 
lower  sign  of  (35b)  and  (39b)  is  negative  as  the  phase  velocity 
is  ,  although  the  magnitude  of  the  group 

velocity  is  much  smaller  than  the  phase  velocity  in  a  range 

ft  ' «  I  .  In  this  range  of  ft'  ,  therefore,  these  waves  ■  are 
highly  dispersive  aa  the  first  inert logravltational  waves  for  roots 
with  the  upper  sign  of  (35b)  and  (39b). 
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6.  Normal  modes 

Arbitrary  waves  in  this  systen  can  be  expressed  as  a  sun 
of  the  waves  of  normal  modes  of  9  »  1  to  6  for  all  order  71  from  0 
to  Infinity,  because  orthogonal  relationships  among  normal  mode 
functions  are  valid  as  proved  In  Appendix  B.  In  order  to  discuss  the 
waves  generated  by  an  Initial  mound  of  the  free  surface  or  by  wind 
stresses  or  atmospheric  pressure,  the  Initial  surface  elevation  or 
forcing  functions  of  wind  stresses  or  atmospheric  pressure  must  be 
expressed  by  a  sum  of  normal  mode  functions. 

The  orthongonal  conditions  for  the  present  system  are 
expressed  as  follows  (see  Appendix  B  for  derivation): 


(44) 


+  VW*  +  ^  $>.»#  +■ 

•  '  *  ap 


in  which  M; 


and 


denote  the  normal 


'i,n9  t 

functions  of  mass  transport  In  x-  and  y-  direction  and  displacement 
(downwards)  of  the  upper  face,  respectively;  the  suffices  1  and  2 
for  j  Indicate  respectively  the  upper  and  lower  layer,  the  suffices 
7l  and  7n  (1  to  6)  and  9  and  f>  (all  positive  Integers)  represent 
respectively  the  mode  and  the  order  of  the  normal  mode  functions  and 
the  asterisk  denotes  the  conjugate  complex.  In  this  equation,  the 
first  and  the  second  three  terms  in  the  integral  sign  represent  the 
sum  of  kinetic  and  potential  energy  of  the  flow  in  the  upper  and 
lower  layer,  respectively  and  tie  last  term  represent  the  effect  of 
interaction  of  the  upper  and  lower  layer. 

The  quantities  Q  niUL')  can  be  computed  by  substituting 
the  normal  mode  functions  for  n=  and  f  =  p  into  the  left  hand  side 
of  equation  (44).  However,  the  values  of  C **(’&')  can  be  estimated 
without  resorting  to  a  complicated  calculation.  As  discussed  In 
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Sections  3  end  4,  the  normal  mode  functions  are  almost  completely 
separated  Into  barotropic  and  barocllnlc  mode  In  a  range  of  0?  much 
less  than  unlnlty  or  In  a  range  of  wave  lengths  from  several  to 
several  thousands  kilometers.  Therefore,  for  the  barotropic  mode. 


the  normal  functions  of  the  lower  layer 
expressed  approximately  by 


Mjl  »  Nx 


can  be 

(45a) 


In  this  case  the  mass  transports  in  the  whole  depth  become 

(45b) 

Hi  '  -h  , 

Substitution  of  (45a)  into  the  left  hand  side  of  (44)  leads  to: 


Ce  *  («) 

+  'N"*' )  j . 

The  integral  equals  to  the  norm  of  the  set  of  eigenfunctions  (or 
normal  mode  functions)  for  a  homogeneous  water  with  the  depth  of 

and  is  determined  by  Reid  (1958)  as  his  equation  (64).  (He 
used  a  notation  instead  of  Cnt  ) 

The  normal  functions  of  the  lower  layer  for  barocllnlc 


mode  are  also  given  approximately  by: 

Mi  ^  -v. 


(47) 


Substitution  of  (47)  into  the  left  hand  side  of  (44)  yields: 


(48) 
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The  Integral  of  tha  first  three  terms  of  the  right  hand  side  is  the 
same  as  the  norm  of  the  normal  functions  for  the  upper  layer  only 
as  treated  by  Reid  (1958) .  The  integral  of  the  last  term  equals 

As  in  case  of  a  homogeneous  water,  the  quantity  : 

p,,„t£)  =  ufcC'oeor  <*»> 

indicates  the  weight  function  of  partitioning  of  the  energy  spectrum 
among  the  six  discrete  modes  (  jf  »  1  to  6).  This  function  also 
satisfies  the  condition  that: 

j;  ps, (&)--/  <*» 

The  weight  functions  for  baroclinic  mode  are  almost  the  same  as  those 
determined  by  Reid  (1958)  for  a  homogeneous  water  with  the 
depth  X).  The  approximate  formula  (48)  shows  that  C*f  2  . 

Therefore,  for  the  baroclinic  mode,  we  have 

?  (51) 

This  relation  indicates  that  energy  spectrum  of  an  initial  surface 
displacement  or  of  forcing  functions  is  distributed  for  the  most 
part  among  three  barotropic  modes,  because  the  weight  functions  for 
baroclinic  mode  are  less  than  $  .  Therefore,  there  is  no  favor¬ 

able  condition  for  exciting  internal  waves  in  a  range  of  wave  lengths 
of  aeveral  thousands  kilometers. 
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7.  Concluding  remarks 

Mathematical  treatment  of  edge  waves  in  a  two-layered 
rotating  sea  reveals  to  us  two  main  features  which  may  be  important 
to  the  application  of  the  theory.  One  is  that  no  instability  occurs 
to  this  system  in  a  range  of  wave  lengths  of  several  to  thousands  * 
kilometers .  In  a  non-rotating  system,  shearing  motion  in  a  verti¬ 
cally  stratified  layer  becomes  unstable  in  a  certain  range  of  wave 
lengths  as  in  Helmholz  waves.  Even  in  a  rotating  system,  a  mathe¬ 
matical  model  of  cyclone  waves  and  a  simplified  theory  of  meander  of 
a  wide  ocean  current  predict  that  instability  occurs  in  a  two- layered 
system  of  geostrophic  currents.  Therefore,  in  this  model  there  is 
no  favorable  range  of  wave  lengths  for  which  some  disturbances  may 
excite  large  internal  waves. 

Another  feature  is  that  energy  of  an  initially  static  mound 
of  water  and  of  disturbing  forces  due  to  atmospheric  pressure  or 
wind  stresses  is  distributed  for  the  most  part  among  three  baro- 
tropic  modes.  This  again  Indicates  that  internal  waves  are  not 
favorably  excited  by  such  disturbances  with  an  order  of  dimension 
corresponding  to  a  range  of  wave  lengths  considered.  This  is  in 
accordance  with  a  general  theory  developed  by  Veronis  and  Stcmmel 
(1958),  who  concluded  that  the  ocean  responds  to  variable  winds 
principally  as  a  homogeneous  body  of  water  and  for  periods  shorter 
than  several  weeks.  Therefore,  possibility  of  favorable  excitation 
of  internal  waves  in  this  system  may  be  found  outside  a  range  of 
wave  lengths  discussed  here.  In  a  range  of  shorter  wave  lengths, 
however,  the  effect  of  Coriolis'  force  is  negligible.  In  a  range 
of  longer  wave  lengths,  the  effect  of  change  of  Coriolis'  coefficient 
may  become  important.  In  either  range  the  present  model  is  not 
adequate. 

Acknowledgement:  Typing  of  the  manuscript  by  Mrs.  Suzy  Bowman  is 
greatly  appreciated. 
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Stoumel  (1953,  1958)  discussed  a  perturbation  problem  for 
a  wide,  geostrophlc  current  Vt In  the  upper  layer  of  two- layered 
ocean.  He  obtained  a  result  that  there  Is  a  narrow  range  of  wave 
numbers  about  ^  =  i/{$  V,  In  which  waves  are  unstable.  In  contrary 
to  the  present  work.  Although  he  assuoed  that  the  perturbation  Is 
uniform  in  x  -  direction  (cross  current  direction),  his  model  may 
be  considered  as  a  modified  use  of  the  present  theory. 

It  is  assumed  that  the  uniform  geostrophlc  current  flows 
In  the  wedgeshaped  upper  layer,  below  which  the  motionless  lower 
layer  extends  infinitely  downwards.  (See  the  attached  figure)  The 

coordinate  system  Is  taken  as  the 
same  as  in  the  text,  but  no  coastal 
boundary  Is  needed  In  this  problem 
Equations  of  motion  and  continuity 
for  the  upper  layer  afe  the  same  as  (1)  and  (2)  of  the  text.  The 
assumption  of  no  motion  In  the  lower  layer  leads  to: 

.  (Al) 

Substituting  (4)  and  (5)  into  (2)  and  using  (Al),  we  have  : 


In  which  07  -  <J~-  -&.V  ,  A,  »  *-  <J7Z  and  the  suffix  1 

for  5  ,  ^  ^  and  y  is  omitted. 

Assuming  that  £  is  independent  on  x  and  that  fl  in  the 

coefficient  of  in  equation  (A2)  is  constant,  Stommel  (1954, 

1958)  derived  a  cubic  equation  of  Q~  by  equating  the  coefficient 

of  ^  to  zero.  However,  if  ft.  is  taken  as  5 1C  throughout  the 

2 

whole  space,  equation  (At)  can  be  written  as: 


5 &  ,  £ 


in  wflitli  UJ  ,  -  t<r-  V 


T 


i-  i1*}  *  - 

t  tAAlA  ' 


0 


(A3) 


The  normal  mode  of  £  which  is  finite  at  x 
responds  to  satisfying 

&U >;*  +  Hi ^  *  -  ran+olW 

This  equation  can  be  transformed  into 

Cc,*-  •+  1^  to,  -  it  =  o 

where  -fc /  =  & V/  ■$  "*  as  defined  in  the  text. 

By  taking  cu,  a  -  (£72  /3  W 


0  and  00  cor* 
(A4) 

(A5) 


,  (A5)  is  changed  into: 


W ’ -t  2.  =■  P  VV 

-^3, 


<A6> 

where  P  -  (*/&'  )  I  ^  (A7) 

when  p<  3  ,  there  are  complex  roots  of  (A6).  However,  as  easily 

seen,  P  £  3  for  all  positive  -&/  .  (The  equality  occurs  for 
k'-2  and  n  -  0) .  Therefore,  in  the  present  model  there  is  no 
range  of  &  in  which  waves  are  unstable. 

When  there  is  a  geostrophlc  current  in  the  lower  layer  as 
in  the  model  of  the  text,  'ft  in  equation  (A5)  should  be  changed  into 
(li-£.)  jk/  ,  in  which  £  »  — "VS/TTi  .  Since  (1+  is  always 

positive  as  discussed  in  the  text,  the  same  reasoning  as  above  is 
valid  and  the  waves  are  always  stable. 

It  is  noted  that  equation  (A5)  is  reduced  to  the  frequency 
equation  (30)  of  Reid's  paper  (1958),  when  aJ  and  are  replaced 
by  —  w($  and  6  $  S  £  " 1  ,  respectively.  Therefore,  mathematical 

features  of  the  roots  of  (A5)  are  the  same  as  those  of  the  barotropic 
mode  discussed  by  him.  Discussions  on  phase  velocity  and  group 
velocity  for  each  of  the  three  barotropic  modes  are  also  applicable 

g 

to  the  barocllnlc  mode,  when  the  velocities  are  interpret' d  as  those 
referred  to  the  coordinates  moving  with  the  flow  in  the  upper  layer 
and  the  wave  number  &  for  the  barotropic  mode  is  replaced  bySttti^. 


Appendix  B 


Equations  (1)  and  (2)  yields  the  equations  of  motion  and  continuity 
for  a  mode  H  and  order  7t  such  as 

.f  -  4 A^„/ »  »  i*1*  2  •  (B1) 


i  M;. v t  rf.,, .  ^  .ml  i_  ' 

-jf—  -nib tj  „e  '  J,  ,,2  . 


(B2) 

(B3) 


If  n  and  1  are  replaced  by  m  and  £,  respectively  and  all 
quantities  are  changed  into  their  complex  conjugate ,  the  following 
equations  are  obtained:  4 

1  '  4  ,  /=  '  •  2  (B4> 

(B5) 
(B6) 


*  7/;.V  j  -  1 ,  2 

d*1i£r  .  /  AN;.*?  *  ' i(T] )=  i,  2 


dz 

where  the  asterisk  denotes  the  complex  conjugate. 

If  equations  (Bl) ,  (B2),  (B3)f  (BA),  (B5)  and  (B6)  are 
multiplied  respectively  by  ,1V, 


and  and  all  the  resultants  are 

summed  up  both  for  j  ■  1  and  2,  the  follotdig  equation  is  obtained 


wf  *  trf  + 


in  which  equation  (14)  is  used. 

When  this  equation  is  integrated  from  0  to  oo  with  respect 
to  x,  the  right  hand  side  of  the  relation  vanishes  because  of 
boundary  conditions  at  x  -  0  and  at  x  •  oo  .  Therefore,  the 
expression  in  brackets  of  the  left  hand  side  of  (B7)  must  vanish 


except  for  n  -  n  end  p  -  1. 

*  StyL  ^  ^  **p  *  N£*p\  +  ^  $>.*«  $».*p  t* 

+■  |j-  (V,  -\/*)fiS>,S  (2H+ |  )qTng Ok<’>f>.(orj<P  •*•  ^T-kC)  -  ty^p^nt  »<C4)  »<r*>»  V 

>2  '  ^flQC-pC^-^tXV^p)  (£$>« 

s  {  0  'h*w\  v*-jaA  &*? 

[c  (it)  fi*m  twd  fi  =  p. 

The  quantity  (^"^can  be  calculated  by  using  the  following  rela¬ 
tions  and  equations  (4),  (5)  and  (14). 


^ f  A(£*_,  t  =o 

£,.r!  -  L-n(2£*) 

),' =  f 


(B9) 


(BIO) 


w  =  n 

W  *  71 


The  result  Is: 


c*p- 


~ — [ ">»' >s.j  £ ;tt,s  rf ')  t its/j  fe  «r»4  A  + 
Z&/*L 


nr 


t  AQhTl^K^  K  fS.s  ft.  f  ia  -5^# 

2  l  <£*  (4  -  <ft«  ;  J + 


(Bil) 


+*'i 
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